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ABSTRACT 

We present ground-based near-infrared (H-band) imaging of the circumstellar disk around the nearby 
classical T Tauri star TW Hydrae. The scattered light image shows a face-on disk with radius 4 arc- 
seconds (corresponding to 225 AU) and a morphology which agrees with recent images from the Hubble 
Space Telescope and the Very Large Array. The best fit power-law for the disk's radial surface bright- 
ness profile obeys the law 7--3-3±0-3. We use our image and published continuum flux densities to derive 
properties of the disk with a simple model of emission from an optically thick, flat disk. The best-fit 
values for disk mass and inner radius are 0.03 Mq and 0.3 AU; best-fit values for temperature, density, 
and grain opacity power law exponents {q, p, and P) are 0.7, 1.3, and 0.9, respectively. These properties 
are similar to those of disks around classical T Tauri stars located in more distant molecular clouds. 
Because of TW Hydrae's nearby location and pole-on orientation, it is a uniquely favorable object for 
future studies of radial disk structure at the classical T Tauri stage. 

Subject headings: stars: individual (TW Hydrae) — circumstellar matter — stars: pre-main-sequence 
— planetary systems: protoplanetary disks — accretion, accretion disks — open 
clusters and associations: individual (TW Hydrae) 



for an optically thick, massive (~ 10^^ ^'^o) circumstellar 
disk. This disk appears to be still actively accreting onto 
the star, as inferred from Ha emission (Muzerolle et al. 
2000). In short, TW Hya shows all the signs of being a 
classical T Tauri star with a surprisingly advanced age. 

Recently, Krist et al. (2000) (hereafter KOO) and Wein- 
berger et al. (1999) reported the detection of a face-on 
(radius « 3.5") disk around TW Hya in HST images ob- 
tained with WFPC2 and NICMOS, respectively. Wilner 
et al. (2000) resolve (radius > 2") thermal emission from 
circumstellar dust around TW Hya with 7 mm VLA ob- 
servations. In this Letter, we present the first ground- 
based, near-infrared image of the circumstellar disk around 
TW Hya, and combine this data with thermal emission 
data to derive disk properties with a simple disk model. 
We find disk properties that are similar to those of disks 
around other cTTs. The nearby TW Hya provides the 
best known opportunity to study the earliest stages of po- 
tentially planet-forming systems. 

2. OBSERVATIONS AND RESULTS 

We observed TW Hya on 7 Apr 2000 (UT) from NASA's 
Infrared Telescope Facility (IRTF) on Mauna Kea. We 
used CoCo, the Cold Coronagraph, to block the light from 
the central star and reduce sidelobe intensity to allow 
imaging of the circumstellar disk; CoCo is used as a front 
end to the IRTF facility camera NSFCAM (Rayner et al. 
1993; Shure et al. 1994). CoCo is a cryogenically cooled 
coronagraph with gaussian apodized focal plane masks and 
a cooled articulatable pupil plane mask (Wang et al. 1994; 
Toomey et al. 1998). We used a mask whose half-power 
(50% transmission) diameter is 1'.'86. A series of eight 
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1. INTRODUCTION 

The formation of the solar system from a nebula of gas 
and dust is strongly supported by observations of disks 
around young stars (Koerner 2001). Observations of these 
protoplanetary disks are largely confined to star-forming 
clouds at distances greater than 140 pc. In these regions, 
the majority of classical T Tauri stars (cTTs) are found 
to have disks with masses in the range 10~^ to 10~^ Mq 
(Beckwith et al. 1990; Andre & Montmerle 1994). Aper- 
ture synthesis mapping of circumstellar gas has established 
that these disks have radii of a few 100 AU (Koerner & 
Sargent 1995; Dutrey et al. 1996), in agreement with di- 
mensions derived from HST images of young disks (O'Dell 
et al. 1993; Padgett et al. 1999). 

The discovery of a nearby lone classical T Tauri star 
by Rucinski & Krautter (1983) affords a unique opportu- 
nity for high spatial resolution studies of a gas-rich disk. 
TW Hydrae (TW Hya, K7Ve, d w 56 pc) is the eponymous 
member of a nearby group of young stars (Kastner et al. 
1997; Webb et al. 1999). The members of the TW Hydrae 
Association have ages of only a few to ten million years 
despite their distance from a molecular cloud (Jensen et 
al. 1998; Webb et al. 1999). TW Hya (the star) stands in 
contrast to its namesakes by virtue of its large Ha equiv- 
alent width (220 A) (Kastner et al. 1997), CO emission 
(Weintraub et al. 1989; Zuckerman et al. 1995; Kastner et 
al. 1997), and molecular hydrogen (Weintraub et al. 2000). 
Furthermore, the considerable 1.1 and 7 mm flux densities 
of 874 mJy and 8 mJy, respectively (Weintraub et al. 1989; 
Wilner et al. 2000) - comparable, after distance scaling, 
to flux densities of cTTs at much larger distances - argue 
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exposures of TW Hya, typically 300 seconds long, were in- 
terspersed with similar ones of nearby sky. Total on-source 
integration time was 2220 seconds. We also obtained six 
20-second exposures of the nearby G star SAO 202496 to 
measure the point spread function (PSF); these observa- 
tions were made immediately after the TW Hya observa- 
tions. SAO 202496, which has a similar H magnitude to 
TW Hya (Amag < 0.5) and a sky separation of 6.4 de- 
grees, is well suited as a PSF comparison star for use in 
analysis of the TW Hya images. We observed both stars 
in H band (Ac =1.62 /xm, AA = 0.28 /xm); seeing was con- 
sistently O'.'S at airmasses 1.7-1.8 for TW Hya and 1.8 for 
SAO 202496. We also observed the photometric flux stan- 
dard SAO 205948 (Ehas et al. 1982) at similar airmasses. 

All exposures were sky-subtracted and flat-fielded. Each 
frame of SAO 202496 was aligned, scaled, and subtracted 
from each TW Hya frame. Alignment and scaling of the 
target and PSF frames were accomplished by minimiz- 
ing the least squares difference (x^) between target and 
PSF frames while varying amplitude scaling and position 
offsets. The least squares difference was calculated in a 
region 4 to 6.5 arcsec from the center of the frame. Image 
shifting in steps of 0.1 pixel was carried out in the Fourier 
domain. The amplitude was varied in step sizes corre- 
sponding to approximately 1% of the best fit value. We 
carried out this same procedure on sets of two independent 
exposures of the PSF star to ensure that PSF— PSF re- 
sults produced a null result and to characterize artifacts 
and residuals. Circular emission 3-4 arcsec in radius was 
detected in all 48 combinations of target— PSF frames with 
a signal to noise ratio of several hundred compared to the 
background RMS and 80 with respect to residuals in the 
PSF-PSF subtraction. Because CoCo rejects 90%-95% 
of the light from the central star, scattered light from the 
circumstellar disk itself dominates from 2" radius out to 
the edge of the disk (see below) . 

Figure la shows our best image of the nebulosity around 
TW Hya. A circular halo of size r=4" (5(7 detection limit) 
is clearly apparent. The peak intensity of the emission is 
6.8±1.1 mJy/arcsec^ at a radial distance of 0'.'94 from the 
star, the smallest radial distance we can see using this focal 
plane mask. The reported uncertainty comprises a conser- 
vative estimate of errors. The largest source of error is 
imperfect centering of the stars behind the coronagraphic 
mask. This effect results in azimuthal variations of ap- 
proximately 15% in the final subtracted image when aver- 
aged over complete quadrants; this is estimated by assum- 
ing that the circumstellar emission is axisymmetric and 
centered on the star. By contrast, the la pixel-to-pixel 
RMS noise in the image background is 10 /iJy/arcsec^, 
corresponding to an error of much less than 1%. (This 
background noise is dominated by noise from the scaled 
PSF star background.) The conservative nature of our 
photometric error estimate is illustrated by subtractions 
of independent observations of the PSF star (Figure lb). 
Here we show the difference between two independent mea- 
surements of SAO 202496; the residuals have a maximum 
above-background emission at just over 1" with an av- 
erage peak surface brightness of 0.1 mJy/arcscc^. These 
images are subject to the same centering errors as the PSF- 
subtracted image of TW Hya, and the residuals seen in the 
PSF-PSF image can represent the error due to imperfect 



centering. The total integrated flux density from the cir- 
cumstellar emission around TW Hya is 21.6±3.5 mJy; by 
comparison, the total integrated fltix density of the resid- 
uals in the PSF-PSF image is around 0.25 mJy. The 
fidelity of the final TW Hya image is bolstered by the 
comparison displayed in Figure Ic. Here we show the 
KOO WFPC2 I band image in grayscale overlain with con- 
tours from our TW Hya image (Figure la). In radial ex- 
tent and orientation, there is a very clear and convincing 
match between the two. 

In Figure 2 we show azimuthally averaged radial pro- 
files of TW Hya and PSF star images before PSF sub- 
traction (thin solid lines) as well as the radial profile of 
the final, PSF-subtracted image of the circumstellar emis- 
sion around TW Hya (thick solid line). Wc also show the 
azimuthally averaged radial profile of the PSF-PSF resid- 
ual (thick dotted line). This figure illustrates a point made 
above: outside of 2", a circumstellar excess is seen in the 
TW Hya image compared to the PSF star image before 
PSF subtraction from either image. This excess extends 
to around 4" in the fully-reduced (PSF-subtracted) data 
(thick solid line in Figure 2), corresponding to an outer ra- 
dius of 225 AU at the distance of TW Hya. A single power- 
law fit to the radial profile follows the law r— 3-3±0-3 (^^jji^jj 
dashed line in Figure 2), although the data does not pre- 
cisely follow a single power law (see Figure 2). Although 
we find scant evidence for the multiple power laws derived 
by KOO, wc do not have the spatial resolution which KOO 
had, and our results are in good agreement with the av- 
erage slope (~-3.5) of the KOO broken power law. Our 
data may show that the radial profile of the disk becomes 
steeper outside of 200 AU, in rough agreement with KOO; 
this level is, however, close to the background level. 

3. MODELING AND DISCUSSION 

The 1.1 mm continuum flux density from TW Hya 

(874±54 mJy) implies a total circumstellar mass 
(gas-Fdust) of order 10"^ Mq (Weintraub et al. 1989). For 
micron-sized grains, this implies a total dust grain cross 
sectional area around 10^ AU^, more than large enough to 
intercept all stellar light if oriented in a spherical shell with 
radius similar to the disk size. Additionally, the fractional 
disk luminosity is ~0.3 (Zuckerman et al. 1995) but the 
ratio of disk reflected light to stellar light is only ~10~^ at 
H band. Both statements imply that the optical depth in 
the disk midplane must be quite high. This precludes the 
use of a model which derives temperature structure un- 
der the assumption of direct thermal equilibrium between 
stellar radiation and individual grains (e.g., Backman et 
al. (1992)). Instead, the temperature structure may reflect 
a variety of effects, including viscous internal heating, re- 
processing of stellar radiation, and super-heating in a sur- 
face layer (e.g., Lynden-Bell & Pringle (1974); Chiang & 
Goldreich (1997)). Without a priori knowledge of the rel- 
ative magnitudes of these influences, we find it simplest to 
freely parameterize the temperature structure in a fit to 
the spectral energy distribution (SED) as in Beckwith et 
al. (1990). 

We derive properties of the TW Hya circumstellar disk 
by using a simple, two-dimensional model which incorpo- 
rates both our imaging data and the thermal continuum 
flux densities shown in Figure 3. Our best disk model is 
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also shown in Figure 3. Following Bcckwith ct al. (1990), 
we assume power-law radial profiles for the temperature 
and surface density, T = Ti{r)-i and E = So(r/ro)-P, 
where r is in AU, Ti is the temperature at 1 AU, and So 
(the surface density at ro) is couched in terms of the disk 
mass, Mdisk, by integrating over the disk surface area. The 
grain emissivity is k = 0.1(z//10^^ Hz)^ cm^ g~^. Ti and 
q are found to be 125 K and 0.7, respectively, by fitting 
points at infrared wavelengths where the disk radiation is 
optically thick; Ti « 125 K is similar to results for other 
cTTs (Beckwith et al. 1990). An inner radius (Rq) of 
0.3 AU is required to fit the 12 /xm point. The face-on 
orientation and outer radius of 225 AU are indicated by 
the image in Figure la. We vary the remaining parameters 
[M^isk, p, and /?) over a grid of appropriate values while 
minimizing the reduced difference between model and 
data. Our derived best-fit values are Mdisk = 0.03 Mq, 
p = 1.3, and /3 = 0.9. 

Typical q values for cTTs (-0.5; Beckwith et al. (1990)) 
can indicate back-warming from an envelope or reprocess- 
ing by a geometrically flared disk (Natta 1993; Chiang & 
Goldreich 1997). Our larger value for q (0.7) is closer to 
the value expected for cither a standard viscous accretion 
or geometrically flat reprocessing disk, perhaps suggesting 
the disk flattening implied for some transitional sources 
(Miyake & Nakagawa 1995). Radial brightness power law 
exponents for geometrically flat to fully flared disks are 
expected to range from -3 to -1.8. Our power-law fit ex- 
ponent is -3.3±0.3, consistent with the average KOO expo- 
nent of ~-3.5. Steeper indices than -3 might result from 
regions of the disk which are shadowed from the star. The 
higher angular resolution KOO observations yield regions 
with power-law exponents between -0.5 to -5.9, perhaps 
suggesting the kind of ripples which can exist for passively 
heated disks (Chiang 2000; Dullemond 2000). 

The inner disk hole [Rq ~ 0.3 AU) implied by our model 
is required by the deficit of disk-derived flux at 12 /xm 
compared to that at 18 /zm. However, more sophisticated 
models of face-on disks can produce a relative decrease 
at 12 /xm without an inner disk hole (Chiang & Goldre- 
ich 1997); our model simply requires an inner disk hole as 
a substitute for more complex radiative transfer calcula- 
tions. Future work includes applying these more sophisti- 
cated models to our TW Hya data. 

We derive a temperature power law index of g = 0.7 
from the disk continuum emission, consistent with the av- 
erage KOO scattered-light disk scale height radial power 
law exponent of ~1.15 (Kenyon & Hartmann 1987). Simi- 
larly, both we and KOO derive (average) density power law 
indices of p = 1.3. Our value for p is not narrowly con- 
strained, and values between 1.1 and 1.5 all provide accept- 
able flts to the continuum data; however, our p is better 
constrained than is typical for SED fits since the principal 



degeneracy is with outer disk radius, which is given here 
by our image. KOO likewise found their density exponent 
not well constrained. A value oi p = 0.75 is the theoreti- 
cal expectation for steady-state accretion (Lynden-Bell & 
Pringle 1974), although Bell et al. (1997) incorporate ef- 
fects such as disk opacity and shadowing and find compli- 
cated radial density distributions with effective average p 
closer to unity. 

Our derived Mdisk of 0.03 Mq and fi of 0.9 are nar- 
rowly constrained by the long wavelength (0.8, 1.1, and 
7 mm) data; Wilner et al. (2000) found Mdisk = 0.03 Mq 
and /3 = 1.0. These derived /3 values differ from /? for 
the ISM (^^2), but agree with that foimd for dust around 
other cTTs (Beckwith & Sargent 1991; Dutrey et al. 1996). 
This effect has been explained as the rcsiilt of aggregation 
of dust grains into larger disk particles (Beckwith et al. 
2000). 

The age derived for TW Hya, 3 to 10 x 10^ yr (Jensen 
et al. 1998; Webb et al. 1999), is somewhat older than the 
ages typically reported for other cTTs with similar disks, 
10^ to 5 X 10*^ yr (Beckwith et al. 1990). It is therefore sur- 
prising that TW Hya possesses youthful disk properties, 
even though it exhibits a stellar age at the upper end of 
this range and disks around other Association members are 
considerably more evolved than the TW Hya disk (Koerner 
et al. 1998; Jayawardhana et al. 1998, 1999a,b; Schneider 
et al. 1999; Koerner et al. 2000). However, because there 
can exist a wide variation in initial disk masses and life- 
times (e.g., Osterloh & Beckwith (1995)) and because ages 
for cTTs in star-forming regions have great uncertainties 
(due to poorly ascertained distances and extinction cor- 
rections), we suggest that disk properties should be taken 
as the key indicator of a disk evolutionary state which is 
only approximately correlated with HR diagram-derived 
stellar ages. In this view, TW Hya is especially useful as 
the closest example of a protoplanetary disk in the classi- 
cal T Tauri phase and with a mass similar to the minimum 
mass solar nebula. Since disk properties at this stage are 
likely to hold great importance for theories about the ori- 
gin of planetary systems, TW Hya's nearby location and 
face-on orientation make it ideal for future high-resolution 
studies of T Tauri disks. 
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Fig. 1. — Images of the TW Hya disk and PSF— PSF residuals. Panel (a) shows our final (1620 second) image. This is a stack of 6 coadds; 
each coadd is a 120 or 300 second TW Hya frame with a 20 second PSF star frame scaled, shifted, and subtracted from it. We have applied a 
linear stretch from (black) to 0.67 (white) mjy/arcsec^ to emphasize the outer radii. The central black circle represents our coronagraphic 
mask and was added in processing; its diameter is equal to the half-power point of the Gaussian focal plane mask. The black circle in the 
lower left is a faint point source in the PSF frames; this point source is at least 10 magnitudes fainter than the PSF star and is therefore 
completely negligible in the PSF subtractions. Panel (b) shows two H band, 20 second PSF images subtracted from each other, and illustrates 
typical residuals in the scaling and subtraction. The linear stretch in this panel is identical to that in panel (a). Panel (c) shows the I band 
WFPC2 image (KOO) with CoCo H band contours (from panel (a)) overlain. The WFPC2 data arc shown in logarithmic stretch (KOO), and 
the CoCo contours correspond to 0.30, 0.65, 1.4, and 3.0 mjy/arcscc'^. The grey and black striped pattern near the edges of the rotated 
KOO image are artifacts of the rotation interpolation scheme we used; the original KOO image has North 18.5 degrees counter-clockwise from 
straight up (J. Krist, pers. comm.). All three panels have North up and East to the left. The field of view is approximately 7'.'04 on a side 
for each panel. 

Fig. 2. — Azimuthally averaged radial profiles for images of TW Hya and the PSF star, both before PSF subtraction, and of the final 
(PSF-subtracted) TW Hya circumstellar disk. The upper and lower thin solid lines show radial profiles of images of TW Hya and the PSF star, 
respectively, before PSF subtraction. The PSF star profile is scaled to be equal to the TW Hya profile at 0''94. The TW Hya profile is 
distinctly different from the PSF profile outside of around 2" (110 AU), showing that we detect the circumstellar disk in the image of TW Hya 
before PSF-subtraction. We also show (thick solid line) the radial profile of the final (PSF-subtiactcd) image of the circumstellar disk (the 
image shown in Figure la). The power law fit to the disk, which follows an r^'^ "^ power law, is shown as the thick dashed line. For comparison, 
the radial profile of the PSF— PSF residual (Figure lb) is shown in the thick dotted line; this residual is quite small compared to the TW Hya 
disk detection. 

Fig. 3. — Spectral energy distribution for TW Hya. The dashed line is a black-body for the star TW Hya (4080 K). The solid line is the 
sum of stellar and disk fluxes, including the best fit disk model found with the five shown parameters varied. All parameters except p are well 
constrained by the multi-wavelength fluxes and near infrared data. Asterisks are near-infrared (JHK) fluxes from Webb et al. (1999); closed 
circles are mid-infrared fluxes from Jayawardhana ct al. (1999b); open circles are IRAS fluxes (IRAS 1986); squares are 0.8 and 1.1 mm data 
from Weintraub et al. (1989); and the triangle is the 7 mm flux density measurement from Wilner et al. (2000). Error bars are small vertical 
lines shown for each data point. 
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